f[x_] := 
        (1/(Exp[((x - 7 10^-19)/(1.69253 10^-23 T))] + 1));
\(p1 = 
      Plot[Evaluate[
          Table[f[x], {T, 5000, 10000, 100}]], 
{x, 
          1\ 10^\(-19\), 
          10\ 10^\(-19\)}, 
PlotStyle -> 
          Table[Hue[0.1  i], {i, 1, 10}], 
Prolog -> 
          AbsoluteThickness[2], 
        AxesLabel -> {"E", "F(E)"}, 
Background -> GrayLevel[0.7]]\)